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SUMMARY 


As a first step toward the computation of the velocity 
distribution along a wing profile of arbitrary shape in a 
compressible fluid, the circulati on- free flow around a sym- 
metrical profile is treated under the assumption of the sim- 
plified dens i ty-soeed relation due to Tchaplygin, Karman, and 
Tsien. The velocity distribution problem is reduced to a non- 
linear integral equation which is solved by a fairly rapidly 
convergent iteration method. Numerical examples are given. 


INTRODUCTION 


The central problem in the two-dimensional theory of a 
potential flow of a perfect fluid around an airfoil profile 

is that of d et ermini ng the pre ss ur e di s t r i but i o n on a profile 

of given share if the speed and direction of the flow at in- 
finity (undisturbed flow) are known. A solution of this 
problem should consist not merely of giving a mathematical 
existence proof but of indicating a method for obtaining nu- 
merical results of reasonable accuracy in a reasonable amount 
of time. 

The difficulty of the problem depends essentially upon 
the pre s cribed speed a t infinity. If this speed does not ex- 

ceed a certain limiting value (depending upon the profile) the 
flow will be everywhere subsonic. For higher values of the 
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speed at infinity the flow becomes partly supersonic (mixed 
or supercritical flow). Finally, it is probable that for 
too high values of the speed at infinity a potential flow 
becomes either mathematically impossible or unstable* The 
case of mixed flow is the more important one, both from the 
practical and theoretical points of view. Nevertheless, it 
seems that the complete solution of the problem of everywhere 
subsonic flows is a necessary prerequisite for a successful 
attack on the problem of transition through the speed of 
sound. (in fact, at present the very existence of mixed 
flows past a profile has not yet been proved.) 

In view of the admitted difficulty of the problem it is 
advisable to develop the mathematical apparatus by consider- 
ing first the simplest possible cases. The most radical sim- 
plification would be, of course, to neglect compressibility 
altogether. Under these assumptions the pressure distribution 
problem has been solved completely. (See references 1 and 2.) 
In the present report the following two simplifying assump- 
tions are made : 

A. Only circulation-free flows around symmetrical pro- 
files are considered* 

B 0 It is assumed that the velocity potential satisfies 
the simplified differential equation resulting from the so- 
called Chaplygin-Karman-Tsi en eauation of state. (Cf. ref- 
erences 3, 4, and 5.) 

Some remarks may be made concerning this second assump- 
tion. In general, the velocity potential cp(x,y) satisfies 
the partial differential equation 


where p is the density of the fluid. Since the density is 
a given function of the speed q 


equation (l) is nonlinear. The function p(q) is determined 
by the pressure-density relation (eauation of state). In an 
isentropic flow the pressure p satisfies the relation 



0 


( 1 ) 



( 2 ) 
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p p - ^ = constant 


( 3 ) 


where Y is the ratio of specific heats for constant pres- 
sure and constant volume. (The standard value of Y is 
1,405.) This implies the density-speed relation 


where a 0 is the speed of sound at a stagnation point and 

p Q the stagnation density. Chaplygin noticed that the 
equation satisfied hy the potential becomes simpler if the 
density -speed relation is taken in the form 


This relation may be obtained formally from (4) by setting 
y = -1. Though this value of Y violates fundamental phys- 
ical laws, it should be observed that only the density-speed 
relation and not the pr e s sur e-d en s i ty relation enters in the 
equation for the potential . 

As a matter of fact, the function (5) behaves qualita- 
tively in the same way as does the function (4) within the 

subsonic range; that is, for 0 ^ q 3 ^ 2a 0 a /(l + Y), and 
for small values of q/a 0 the function (5) gives a good 
numerical approximation to (4). 

Von Karman and Tsien justify the use of the value Y = -1 
by the remark that it is possible to determine such values 
of the constants A and 3 that the pressure-density rela- 
tion 


will give a good approximation to the relation (3) for val- 
ues of p and p close to some preassigned values, say to 
the values of p and p for the undisturbed flow. This 
remark is of interest as far as computations of the pressure 
distribution are concerned. It in no way affects the velocity 



y/Cv- 1) 


(4) 


1 



(5) 


p = A/p + B 
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distribution for, as it was already noticed, the differential 
equation for the potential function depends only upon the 
density-speed relation, and the preceding pressure-density 
relation leads to the same equation (5) no matter what values 
are assigned to A and B. 

It should he emphasised, however, that the primary pur- 
pose of this report is not to facilitate the use of the ap- 
plication of the approximate relation (5) hut rather to de- 
velop methods which could he extended to the case of the ac- 
tual dens i ty -speed relation. 

In the following, use will he made of certain results 
contained in a previous report. (See reference 6.) 

This investigation, conducted at Brown University, was 
sponsored hy and conducted with the financial assistance of 
the National Advisory Committee for Aeronautics. 

The author largely profited from several instructive 
discussions he had with Professor S. E. Warschawski. He 
also is indebted to Mr. Charles Saltzer for competent assist- 
ance. 


SYMBOLS 


A(tu) auxiliary function defined hy equation (35) 
a local speed of sound 

a 0 speed of sound at a stagnation point 

B(uu) auxiliary function defined hy equation (35) 

C, Cj positive constants 

dS 3 non-Eucl idean length element defined hy equation (22) 

3(P) domain exterior to the profile P 

F integral transformation defined in section 5 

f (cu ) function defining the mapping of the circle into the 
profile P 

f n (ui)) nth approximation to the function f('Jj) 
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G complex potential of a compressible flow 

g(a) function inverse to f(ou) 

Sik coefficients of the metric (22) 

h(uu) function defined by equation (44) 

Im( ) imaginary part of ( ) 

P profile surrounded by a compressible flow 

p pressure 

q local speed 

goo speed of the undisturbed flow 

q(<?) value of q at a boundary point 

q* distorted speed 

a* distorted speed of the undisturbed flow 

00 

a*(a) value of q* at a boundary point 


M local Mach number 

stream Mach number 

R radius of the circle in the £-plane 

Re( ) real part of ( ) 
s arc length measured along P 

S length of the curve P 

t parameter occurring in section 8 

u,v components of the velocity 

w * distorted velocity 

x,y Cartesian coordinate in the z-plane 

z complex variable 
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Z(s) coordinate of the profile as a function of the arc 
length 

Z L leading edge 

Zj trailing edge 

a angle at the trailing end 

y exponent in the adiabatic relation 

e thickness parameter of a symmetri cal Joukowski profile 

e i constants occurring in section 8 

0(a) slope of the profile P 

6 angle between the velocity vector and the x-axis 

6 value of 0 on the boundary 

£ auxiliary complex variable 

A square of the distorted speed of the undisturbed flow 

A function defined by equation (54) 

h Cartesian coordinates in the £-plane 
p density 

p 0 stagnation density 

a dimensionless length parameter along the profile P 

cp velocity potential 

qj value of cp at the boundary 

X auxiliary analytic function defined by equation (34) 

stream function 

u> argument of a point on the circle |£| = R 
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ANALY SIS 

1. The Boundary Value Problem 


Consider a symmetrical profile P in the plane of the 
complex variable z = x + iy. It will be assumed that P 
is a smooth curve, except, perhaps, for a sharp angle at the 
trailing edge z <p , that the x-axis is parallel to the axis 

of symmetry of the profile and that the profile is given by 
an equation of the form 

z = 2(s), 0 < s < S (6) 

where s is the arc length on the curve P measured in the 
counterclockwise direction from the point ztp. Then S is 
the total 1 engt h of the profile and 

z L = Z (S/S) 

is the leading edge. It will be convenient to introduce the 
dimensionless parameter 


ct = 2tts/S (7) 

The function 

©(a) = arg Z'(aS/2TT) (8) 

whe r e 

Z 1 ( s ) = dZ/ds 

depends only upon the shape but not upon the size or posi- 
tion of P. Note that by virtue of the foregoing assumptions 

0(0) = tt - a/2, 0(tt) = 3tt / 2 , 0(2tt) = 2tt + a/2 

where a is the angle at the trailing edge, 0 "5 a 5 tt, 
and 

0(2tt - a) = 3tt - ©(ct), 0 ^ a $ tt (9) 

The equation of the curve P may be written in the form 
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z 


Z ( a S / 2tt ) 



da + 


Z j 


( 10 ) 


Now let aoy(x*y) be the potential of a circulation- 
free flow of a compressible fluid past the profile P; that 
is, a function such that 


u = a o 



v = 


a 0 


by 

dy 


are the components of the velocity in the x- and y-directions, 
respectively, a 0 "being the speed of sound at a stagnation 

point. -he function y(y,y) is defined and one-valued in 
the domain E(P) exterior to P and satisfies the “boundary 
c ondi t i on 

« 0 on P (11) 

dn 

as well as the condition 

dc P a co dc P , v 

— — > — » — — > 0 as z — > oo (12) 

dx a 0 dy 


Here d/£n denotes differentiation in the direction normal 
to P, ar.d q ro is the speed of the flow far away from the 

profile (undisturbed flow). 


The conjugate complex velocity is given by 


< - ie 
w = u - iv = qe 


(13) 


where q is the speed and 8 the angl e between the velocity 
vector and the x-axis. The function 9 satisfies the condi- 
tion 


0 = 0 - tt on the upper bank of P 

0 = © - 2tr on the lower bank of P 


( 14 ) 


and 


0 — > 0 


as 


z 


00 
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Under the assumption of the approximate density-speed 
relation (cf. Introduction) 


P = 


Jl + q 2 /a 3 


the equation of continuity 


JL 

dx 


(pu) + ~ (pv) = 0 
dy 


takes the form 

ii. 2 dtp ijL + r i + (i£f\ afsE. = 0 (15) 

l ^-dy/ J dx 3 dx dy dxdy l_ '•dx/ J dy 3 

This is the classical eouation of a minimal surface. 

The determination of the flow around a given profile P 
requires the integration of the differential eouation (15) 
\inder the boundary conditions ( 1 1 ) and (12). In the case of 
an incompressible flow the corresponding boundary value prob- 
lem can be reduced to the problem of mapping the domain E(P) 
conformally into a domain exterior to a circle. A similar 
mapping will be defined presently for the flow considered 
here . 


2, Mapping of the Profile into a Circle 

The stream function of the flow \|/(x,y) is defined by 
the equations 

dqp ^ p 0 dqp 
dx p dy 

dgp _ P 0 d\|/ 

dy p dx 


This function is constant along any streamline and can be 
normalized so that 
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\|f = 0 on P 

The complex potential G(z) is defined by 

G(z) = cp(x , y ) + i\|/(x , y ) (16) 

Let the potential cp be normalized so that 

cp L = 

where cpj, denotes the value of cp at z = z^ and <Pip the 
value of cp at z = s^. This can always be achieved by add- 
ing a constant to q? , The function 

G = G(z) (17) 

maps the domain E(P) into the domain in the G-plane ex- 
terior to the slit 

V * 0, -cp T - cp ^ cp T (18) 


This latter domain is now mapped conformally into the domain 
l£l - R in the plane of the complex variable £ = £ + ih 
by means of the relation 


Equations (l?) 


of the domain 


are taken into 



and (19) define a transformation 

t = l(x , y ) , h s p (x,y ) 

E(P) into the domain |£| - R. 

Z = Jj t Z ~ 2 |j t 2 a co 

the points 


(19) 


( 20 ) 

The points 


£ = -R , £ = R , £ = oo 


respectively. 


If R is chosen as 
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cp^ao 

2a_ 

* oo 


the mapping (20) satisfies the conditions 


( 21 ) 


di 

dh 



1 


as 


( 22 ) 


By virtue of the foregoing mapping there exists a one- 
to-one correspondence between the points of the profile P 
and those of the circle |£| = R. This correspondence can "be 
described by means of a function 


a = f(a>) 

such that the point Z [f (tu) S/ 3tt] corresponds to the point 

Re lvU . Plainly f(uu) is an increasing function satisfying 
the conditions 

f(O) = 0, f (tt ) = tt, f ( 2tt ) = 2tt (23) 

as well as the symmetry condition 

f (2n - oi) = 2tt - f (uu) , 0 < tu < tt (24) 

In the following sections it will be shown that the knowledge 
of the function f(u>) implies the knowledge of the velocity 
distribution along P. 

Remark "• In the case of an incompressible flow the mapping 
just constructed is exactly the standard conformal mapping of 
the profile into a circle. In the case considered here the 
mapping (20) is conformal with respect to the Riemann metric 
(dS) defined in E(P) by means of the formulas 

dS<! = &iidx 2 + 2g 13 dx dy + gs^dy* (25) 

where 


S 1 1 = 


a 0 3 + q 3 cos 3 0 


a 0 + q 


3 
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Sis ~ 


sin 0 cos 0 

„ 3 , 2 

a o + <1 


Sss ~ 


2 . 2 .2a 

a 0 + q 3 i n 0 

2 . 2 

a Q + q 


The proof of this assertion follows immediately from the re- 
sults of a previous report. (See reference 6.) 


3. Velocity Distribution Expressed in Terms 
of the Function f(uu) 

At a point £ = Re 1 equation ( 1 9 ) takes the form 


or, by (21) 


Cp = Cpip cos u> 


cp 



cos 


UJ 


(26) 


Now let q(cr) denote the value of q at a point Z[aS/2TTj 
and <p(a) the value of cp at this point. Furthermore let 

u> = g(o) (27) 

be the function inverse to f(uu). By (26) 


^ o 

cp(cr) = 2R — — cos g(a) (28) 

a o 


On the other hand, on 

the profile P 





dip (a) 

2na n 



q(a) = 


ds 

S IH ^ 

(a)l 

(29) 

so that 






q(o) = 

4ttR i 

s a -» 1 

s in g (a ) | g 1 

(a) 

(30) 
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This formula shows that the function g(or), and therefore 
also the function f(u)) determines the velocity distribution 
along P but for a constant factor. 

A formula permitting a complete determination of the 
velocity distribution can be derived by introducing the so- 
called distorted velocity w * defined by 

w* = q*e“ 10 (31) 

where a* is the distorted speed given by 


q* = - (33) 

a o 1 + /l + q a /a 0 3 

Note that q* always satisfies the inequality 0 < q* < 1. 

It has been shown (see, for instance, reference 6) that the 
complex potential G is an analytic function of the variable 
w *. Therefore w * is an analytic function of G and hence 
also of the c ompl ex variable £ . The function w * does not 
vanish, except at the points £ = -H and £ = R. The imag- 
inary part of the logarithm of w* is -9. Along the circle 
l£| = R the fxinction -9 may be regarded as a function of 
the real variable (JU, £ = Re^. This function possesses 
jumps of the magnitude a and tt at uu = 0 and U) = tt , 
respectively. It follows from known theorems of function 
theory that at £ = -R the function w* vanishes a 3 (£ + R) 

and at £ = R as (£ - R) 0 ^ 11 , respectively. Furthermore 


w*(co) = q^> 0 (33) 

Hence the function 

x(U - »*<(, * R)- 1 (t - R)' 0 ^ t 1 ^ (34) 

is regular for |£| > R, continuous for f£| = R, and no- 
where equal to either zero or infinity. 

Therefore log x(£) is a one-valued analytic function 
which is continu ou s on the circle |£| = R and regular at 

£ = oc. Set 
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log x(Re^ U ) = A (id) + iB(a’) 


(35) 


A(u>) may be expressed in terms of 3 ((d) "by means of the 
well-known formula 


TT 


A (id) = - 


H 

~ I <! 3 (id +t ) 
n K 0 


B(oi - t ) 1 cot f dt + A 0 

J 2 


(36) 


where A^, is the value of log !Xi at infinity. (See, for 
instance, reference 7, p. 243.) Now "by (33) and (34) 


A» = log q, 


A (id) = log q* [ f (tu ) ] - log |l + e lUJ ! il - e 


a 

iCDI 17 


(37) 

(38) 


and 


B (id) = — 6 [f ( u)) ] — a rg(l + e* °) - ~ arg(e iUJ - 1 ) + ^1 + ^ ^(D (39) 

Here 6(a) and q*(a) denote the values of 6 and 
respectively, at the point Z(aS/2rr) of P. Noting 


q . 

that 


1 + e lU) l = 2 


cos 


and that 


arg (l + e ia ) =• 


( iuj i \ us . tt 

arg(e - 1) = - + - 

Ct (J 


UJ I 
2 i 

1 - e 

r 


tjy 

2 

f or 

i*- 

for 


o . OJ 

2 s in — 


for 0 < id < tt 


as well as that by (14) 

0[ f (id) ] - TT 
0[f(u>)3 - 2 tt 


6[f (uj)J 


r 0[ 
\ ©r 


for 0 < (JU < TT 

for tt<ud< 2tt 


(40) 
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equations (38) and (39) can "be written in the form 


log q* 


f(uj) 


= A(uj) + log < 2 


i + * 
TT 


COS — 


sin 



B(uj) = -0 


f (uj) 


a + tt 

+ 2tt lU + 


(•-t) 


From (41), (42), and (36), it follows that 


f (w) 


1+ Si 




TT 

UJ 


U) 

= C 2 

cos g 


sin ^ 


h (uj) 


where 


TT 





0 


-© 


f (uj - t ) 


f (U) + t ) 

Since hy (32) 

— = - (— - q*) 
q 2 Vq* J 

it follows from (43) that 


— - t > cot ~ dt 
TT 2 


a. 


f (uj) 


a + a, 
q * tt 

«> 2 


UJ 


ID 

'COS — 


sin — 

1 2 


2 


-h (uj) 


3(1 + fr) 3 (ju 
-A2 cos ~2 


where 


UJ 

sm g 


£a 

n h(ou) 
e 


A = 


U*> 

oo 


2 q. 


1+/1 + 


ao 2 


(41) 

(42) 

(43) 

(44) 

(45) 


(46) 


2 


(47) 


16 


NACA TN No. 1006 


This is the desired expression of q(o) in terms of the 
function f(ui). The parameter X may he used instead of 
ci-a>/ a o to determine the conditions at infinity. This param- 
eter can he easily expressed in terci3 of tne stream Mach num- 
ber (of. reference 6): 


X z. 



1 + J 1 - 


( 49 ) 


The fact that tho velocity distribution can be expressed 
in terms of the function f(uu) in two different ways permits 
the derivation of an integral equation for the function f(o)). 


4. The Integral Equation for the F’anction f(uj) 
Equations (28) and (29) may be written in the form 


cp 


q If (u) ) 
L 


f(w)] . 

q„. 

23. — — cos 

0 ) 

J 

R o 


3ua c 

d$[f (UJ ) ] 

1 

S 

d uu 

F'lwT 


Combining these two equations yields the relation 

^ ( l C o 


, / , 2tt a. 

f (oj ) = 5 


S q[f(o))] a 0 


| s in U)| 


(49) 


Now substitute in (49) the value of a 0 /q given by (46). 
Then 


f ' (ui) = C 
where 


uu 

sin 


.a 

'u c 

[e 






- h(u, >-A2 


COS 


ii i 


UD 

,in 2 


(50) 


uR q. 


S a 0 q 


* 

00 


and h(o>) and X are given by (44) and (47), respectively. 
Integrating 
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f (UJ ) 


= C 


UJ 


sin 


0) ' 


■ X2 


TT ' 


r e(n-Si) 

n e 

'v. 


COS 


ur 


sin 


(JL ) 1 


TT 




Setting a) = 2 tt here it follows from (24) that 
2 tt 


,2TT 


= C 


sin 


(JU 


a - 


a, 


2a 




(JU 

cos 


sin 


(JU' 


" e b<u ' l) jda J ' 


so that finally 
f (oj) 


= 2tt- 


UJ m 

J 'sm,,~ 
c i <*• 


ii-fi 
* I TT 


<> + S)i 


~ h (w )_ A2 


1 

2 

, . 1 

, ~ UJ 

| C0S '2~ 


sinyi- 

2 


’ ,a . h(u ">jdu,' 


/ 


2TT 


s i n -- 


L'J 


i|'^| e - h ta , ). X2 2 ( 1+ tr' ) 


(51) 


W f 

COS- 


sin 


(JU 1 


2 a 


"e h( “' ,, ju U ' 


Sinco h (uO is given "by formula (44) this is a nonlinear 
integral eauation for the unknown function f(ui). 


5. Solution of the Integral Equation 
The integral equation (51) can he written in the form 

f (ou) = F jou, f (u) f ) j 

where p|oj, f (w f )j> denotes the right-hand side of (51)!. The oper- 

ation F is a functional transformation which takes a con- 
tinuously differentiable function f(u) 1 ) satisfying the con- 
ditions 

f(0) = 0 f ( 2 tt ) = 2tt 


into a function satisfying the same end-point conditions* 
Therefore the solution of (51) can be attempted by the iter- 
ation method. Choose some function f 0 (uj) satisfying the 
preceding conditions and compute successively 
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f a (u>) => r 



fi (vu) 


} 


' n+ 1 


(oj) = F-{uj, 



If the sequence 


f 0 (uu) , f 1 (to) , ... f n (to) , ... 

converges toward a function f(uj) and lim F(f n ) = ( f ) . this 

function f satisfies the integral equation. 

From the purely mathematical point of view it would he 
necessary to supplement the preceding consideration hy prov- 
ing that under suitable assumption: (i) the integral equa- 
tion possesses a solution, (ii) this solution can he obtained 
hy iterations, and (iii) this solution is an increasing 
function. It is hoped that such proofs will he presented at 
come later date. At present it may suffice to state that 
the statements (i) to (iii) seem to he verified in the cases 
for which the computations have been carried out. The ex- 
istence of an increasing function satisfying the integral 
equation seems quite obvious from physical reasons. As for 
the convergence of the method, reference is made to the fact 
that the method described here is rather similar to 
Theodorsen's method of conformal mapping (references 1 and 2) 
for which a rigorous convergence proof has been found (ref- 
erence 8 ) , 

It might he noted that the desired solution f(cu) must 
satisfy the symmetry condition 

f ( 2rr - uu) = 2tt - f(uu), f’(2n - o>) = f'(uj) (52) 

If the function f Q ( UJ ) satisfies this condition, so will 
all successive approximations f n (uu). It will therefore he 
sufficient to compute f n (o>) only in the interval 0< UJ< tt. 
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The only nontrivial step in computing the functions 
f n (uj) consists in evaluating the integral 


TT 


n+i 
wher e 


2tt , 



h 4 - ( UJ ) = ~ / <A (a' + t ) - A , (u) - t ) } cot p- &t 


n+i 


A n + 1 (oj) = 0 


n + i 


f n (a)) 


a + tt 
2tt 


0 ) 


(53) 


(54) 


(Cf. equations (5l) and. (44). ) It should he noted that this 
is a proper Rienann integral. In fact, the value of the in- 
tegrand at t = 0 is 


1 iia 

t-»o 



(oj + t) - A n+ 1 
2t 


(u) 



cot 


t, 

2 


J 


= 4Ai +1 (tu) = 4 j©' Cf n (cu)] f£(u>) - 

By using this information, the integral (53) can he evaluated 
numerically, say hy the trapezoidal rule. 


After f On ) has been computed with sufficient accuracy, 
the velocity distribution is computed hy means of formula 
(46). 


6 . Choice of the Function f 0 (uu) 

The rapidity of convergence of the iteration method for 
solving equation (51) will depend upon the choice of the 
function f 0 ( ai ). the Oth approximation. In order to reduce 

the computational work, this function should always satisfy 
condition (52). 

A few methods of choosing the function f 0 (uu) are 
listed, in the order of preference? 
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(а) Choose for f Q (u)) the solution of equation ( 51 ) 

for a value X* as close as possible to the value of X for 
which the equation is to be solved. 

(p) Choose for f 0 (uo) the solution of equation (51) for 

the desired value of X or for a value X' close to the de- 
sired value, and for a profile P 1 different from but close 
to P. 

( Y ) Choose for f 0 (tu) the function resulting from the 
conformal mapping of the profile P onto a circle; that is, 
the solution of (51) for X = 0, 

(б) Choose for f 0 (uu) a function approximating the 

function resulting from the conformal mapping of the profile 
P onto a circle. For thin profiles such a function is given 
by 

f 0 (uj) = ; (1 - cob «)), 0 ^ uj < it 

* (55) 

f 0 ( 2tt - u>) = 2tt - f 0 (u>) 

Note that (Y) is a special case of (a) (set X* = 0) and (fi) 
a special case of ( P ) (set X 1 = 0 and choose P as a 
straight segment)* 


7. Velocity Distribution at Points Not on the Profile 

It remains to show how the knowledge of the function 
f ( uj) permits the computation of the velocity distribution 
at points not on the profile. This is done by means of the 
following theorem which also shows that solution of the in- 
tegral equation (51) actually yields a solution of the bound- 
ary value problem stated in section 1* 

Note first that from the way the integral equation has 
been set up it follows that there exists an analytic function 
w*( £) regular for | £ | > 1 and such that 


f (w) 


-16 [f (<o) 3 

o 


w * ( e iuj ) = q* 


(56) 
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where q* is given by (43) and. (44) and 6 by (40). If 
f(e luJ ) is known, w*(£) can be computed, say by Cauchy r s 
f ormula J 


w 


(U = - 


2tt , , 


n 


2TT 


w * ( e ^ 1-0 do) 


1UJ 


- i 


(57) 


Now the following theorem holds: 

I'et f (oi) be an increasing function satisfying (51) . 

Set 


z 


= x + iy = 0. 



t 

**U) 



(58) 


whe r e C 1 is a rea l constant and th e bar denotes the con - 

jugate complex quantity. The transformation 


x = x(£,ri), y = y(£,h) (59) 

of the z -plane into the £ - plane defined by (58) (for | £|>l) 
is one-to-one . It takes the domain | £ [ > 1 into the domain 
E(F) exterior to the profile P. The function 


qp = SC^Re (t, + j) (60) 

considered as a function of x and y is the desired poten - 
tial of the comp re ssible flow around P ; that is. it satis - 
1 e s the differential equation (15), the boundary condition 
(ll) and the condition (12) . 

The proof of this theorem will be found in the appendix. 

After cp is found, the velocity components u and v 
can be determined by differentiation. But it is also true 
that w* considered as a function of z is the distorted 
velocity (cf. sec. 3) and therefore 


22 


NACA TN No. 1006 


. 2 |w *1 w * 

u - iv = a n — 

1 - |w*i 2 I w ^ | 

The proof of this last statement is left to the reader. 


8. Examples 

As an illustration of the method, velocity distributions 
have been computed for a circle and for a symmetrical Joukowski 
profile ’vith e = 0,15 (e being the usual parameter determin- 
ing the thickness). The following values of A have been used 

A = 0.045 for the circle 
A = 0,157 for the Joukowski profile 

These correspond to the following values of the stream Mach 
numb sr 


= 0,406 for the circle 
= 0.685 for the Joukowski profile 

These values of are known to be close to the critical 

values, (The critical stream Mach number is the stream Mach 
number for which the maximum local Mach number is equal to 
unity . ) 

In the case of the circle 

©(a) = a + ©’(a) = 1, a = tt 


It is natural to set 


f 0 (cu) = ui 

This corresponds to case (Y ) of section 6. The first approx- 
imation is easily computed in closed form and is equal to 

f i (uj ) = uo + — — sin 2 U) 

1 - 3A 


In the case of the Joukowski profile 
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a = 0 

The function © ( s ) and ©^s) are given by the parametric 
f ormul as 


2c j (l + e ) f , . . , , ^ 

= n * — ^ e 3 tanh (e 3 cos t ) - c 2 tanh ( € 3 cob t)> 


4 cot t 


0=tt + t - tan 


-l 


(cot t) 2 + — , . 2€ 


/ ,2 l + 2 e 

(cot t ) + 


-|2 


- 4 ( c ot t ) 


3 €* 


t) ‘ + 6(sln _ 

ds l c 3 ( 1 + t ) a 

where the parameter t ranges from t = 0 to t » H- and 

2 

€ 2 > and C 3 are constants determined by 


c 


i 


21£ 


2(1 + € ) 


tanh" 1 c 3 


c 2 tanh” 1 € a 


€ 


2 



e 


3 



+ €)(! + 3€) 
1 + 2e 


ihe proof of these formulas will be found in reference 9. 

The function s, 0, d©/ds are tabulated in table I. The ap- 
proximation of order 0 has again been chosen according to 
case (y) of section 6. In the case of a circle the func- 
tion f(u)) must satisfy the symmetry relation 


f (tT - 0)) = TT - f (OD ) 
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It Is therefore sufficient to consider this function in the 
interval 0 < uu < tt/2. Accordingly the functions f n (u)) 

have teen computed for u» = 0°, 10°, 90°, In the case 

of the Joukowslci profile, the functions f n (ui) have been 
computed for <JU = 0°, 10°, 180°. The convergence of 

the successive approximations is seen from tables II. 

The resulting velocity distributions are given in tables 
III and plotted in figures 1 and 2, The argument 5 i3 the 
argument of a point on the circle into which the profile is 
mapped conformally. The results obtained have been compared 
with those arising from the Karman-Tsi en velocity correction 
formula 

X = f±.\ -L (61) 

1 - 

00 x 

/ q v 

where f — j is the value of the ratio at local velocity to 

\q 

velocity at infinity for an incompressible fluid. To use 
this formula amounts to replacing the function f ( uu ) .by the 
function arising in conformal mapping of the profile into 
the circle. In the case of the Joukowski profile Kaplan's 
results (reference 10) obtained by a modified Poggi method 
are also given for the sake of comparison. 

It will be noticed that the present method (which con- 
sists of an actual solution of the boundary value problem 
for the esse 7 = -1 ) gives a greater compressibility effect 
than the one predicted by the approximate methods mentioned, 
(To evaluate this remark correctly, note that Von Karman ex- 
pressed the opinion that in the case when the assumption 
7 = “1 is applied to air and formula (61) is used, the error 
committed in using this formula seems to counteract the error 
committed in using the incorrect pressure-density relation.) 


CONCLUDING REMARKS 


It has been shown that under the assumption of the lin- 
earized pressure-volume relation and of a symmetrical flow 
the velocity distribution of the compressible flow past a 
wing section of arbitrary given shape can be determined rig- 
orously by a method which requires not considerably more 
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computational labor than the case of an incompressible flow. 
This is, of co\irse, only the first step toward the complete 
solution of the velocity distribution problem. The next 
step should consist of extending the present method (a) to 
the case of the actual adiabatic pressure-density relation, 
(b) to the case of a circulatory flow around a not necessar- 
ily symmetrical obstacle. 

Remark : After this paper was completed the author 

learned about a paper by Slioskin (reference ll), in which 
the same problem is reduced to an integro-dif f erential equa- 
tion, different from the one derived in this paper. 


Brown University, 

Providence, R. I., Kay 1945. 
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appendix 

This appendix contains the proof of the theorem stated 
in section 7, 

The mapping properties of the function (58) follow im- 
mediately from the following three statements. 

(a) The function (58) takes the point l - "° into the 
point z - co , 

(b) The function (58) maps the circle = 1 in a one- 

to-one manner into the profile P. 

(c) The Jacobian 

j _ d »y ) 
d (i ) 

is positive for all values of i and D , | 2 + P 3 > 1. 

To verify ( a ) . observe that as £ — 5 00 , w* approaches 
the value q ¥ = v/a > 0. 

To verify (b), note that the integrals in (58) are inde- 
pendent of the path since the integrands are analytic func- 
tions of £. In order to obtain the image of |£| = 1 the 
integration may be performed along the circle. But for 

£ = e luJ 
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2 sin u) d ut) 


whereas w*(e iaj ) is given "by (56). Hence the equation of 
the curve into which |£| =1 is taken may He written in 
the form 


= 2 C 1 


r 


uu 


i 


•q*[f (UJ)] 


- q. a *‘[f*(oj)3 sin uue 


i £rr+ 9 [f (uu) 


duu 


■jet o = f(u)) and note that the integral equation (51) im- 
plies that f ( 2tt ) = 2rr . By virtue of (40), (45), and (49) 

the preceding equation may he written in the following forms 


z 


C 3 



f ' (uu) 


i©[f (uj) ] 
e dui, 


0 < uu < 2tt 


(A1 ) 



0 


01: ct < 2tt 


( A2 ) 


where C 2 is a new positive constant. Equation (Al) shows 
that l£l = 1 is taken into P (of. equation (lO) and note 
that P is determined hut for a scale factor). Equation 
(A2) shows that the mapping of the circle into P is one- 
to-one, for by hypothesis, 


f'(w) > 0 (A3) 

To verify (c), observe that it follows from (A3). (43) 

and (50) that |w*(e^ aj )| < 1, Since the maximum of the mod- 
ulus of an analytic function is attained on the boundary, 


I w*( 0 | < 1 for j £ J - 1 ( A4 ) 


Now the Jacohian is equal to 
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j . i« (51 . i„ A ( i* + |.iv|e . 5£\) 

\bi br\J \ Vd£ d J / vd J b{J j 


- a »’ u ( 1 -T?>'"( 1 -f)j & ( l -ir>“-( 1 -Ti 



= Cl 1 - 


so that "by (A4) 


v; ■ 


"*D 


J > 0 for |t ! > 1 


Next, equations (58) and (60) may be rewritten in the 


form 


wher e 


= HeXiU), y = Rex a (£), z = Hex. (S') 


Xl(0 = 0 l /{^tlT- w * a) } 0 -F) i£ 


(A5) 


a 


<■ 

/>£ 


*•«> ■ -».y {^Tu + a,<£) }( 1 -F) d£ 

R 


(U = 20 a (t + 1 ) 


Since 


IS is is 

Xj + Xg + X 3 = 0 


(A5) is a We ier s t ras s i an parametric representation of a mini- 
mal surface. In other words, cp considered as a function of 
x and y satisfies equation (15). 
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A Bimple computation shows that a line element normal 
to the circle | £ I = 1 is taken hy the mapping (58) into a 
line-element normal to the profile P. Since the normal 
derivative of cp in the £-plane vanishes so does the normal 
derivative of cp in the z-plane. Thus cp considered as a 
function of x and y satisfies the boundary condition (ll). 

Since 

dz dx dy b z dz 

5 i ~ b i bi ‘ d £ b~t 



and w*(oc) = 


q* . 

00 


dx 


it follows from (58) that as 
1 - q* 8 


q* 


OD 


by 

n 


0 


so that 



dy 

dn 




1 + 


q 


co 



1 q* w 

~ T » “ 

dx Cj 1 - q* s by 


— o — q » 

dx * dy C 1 1 + q* 3 

* CD 


Now, as t, 0 


and theref ore 


dcp 

— ~^2C 1 , 
b i 


dcp 

b r\ 


5<p _ dcp d|_ + dcp dh _ 3q So _ qro 

dx d t dx db dx 1 - a a» 

dcp _ d£ dj_ + dt£ dh ^ 
dy d | by dh by 

as z -$► c». Thus condition (12) is also verified. 
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Functions entering into the computation of the velocity 
distribution along a symmetrical Joukowski profile with 

6 = .15 


t 

s 

@ 

d©/d*| 

T~t 

s 

® 

d©/dS 

0.00 

0.0000 

3.1416 


0.80 

3.0267 

3.8755 


.02 

.0667 

3.0834 

-.4076 

.82 

3.0320 

3.8980 


.04 

.2548 

3.0353 

-.1731 

.84 

3.0371 

3.9204 


.06 

.5191 

3.0034 

-.0804 ! 

.86 

3.0419 

3.9428 


.08 

.8148 

2.9888 

-.0220 j 

j 

.88 

3.0465 

3.9650 


.10 

1.1062 

2.9892 

+.0238 1 

.90 

3.0510 

3.9873 


.12 

1.3723 

3.0009 

+.0647 

.92 

3.0551 

4.0094 


.14 

1.6066 

3.0206 

+.1040 

.94 

3.0591 

4.0315 

1 

.16 

1.8071 

3.0453 

+.1437 * 

.96 

3.0629 

4.0535 

+5.8608 1 

.18 

1.9766 

3.0730 

+.1850 | 

.98 

3.0666 

4.0755 

+6.0844 j 

.20 

2.1193 

3.1024 

+.2287 i 

1.00 

3.0701 

4.0975 

+6.3066 1 

.22 

2.2395 

3.1326 

+.2757 | 

1.02 

3.0735 

4.1194 

+6.5277 •• 

.24 

2.3407 

3.1630 

+ .3264 

1.04 

3.0768 

4.1412 

+6.7464 .| 

.26 

2.4265 

3.1933 

+.3813 

1.06 

3.0800 

4.1630 

+6.9630 ' 

.28 

2.4995 

3.2232 

+.4408 

1.08 

3.0831 

4.1848 


.30 

2.5620 

3. 2526 

+.5052 

1.10 

3.0861 

4.2065 


.32 

2.6158 

3.2816 

+.5750 

1.12 

3.0890 

4.2282 


.34 

2.6623 

3.3101 


1.14 

3.0918 

4.2499 


.36 

2.7029 

3.3381 


1.16 

3.0945 

4.2715 


.38 

2.7385 

3.3655 

+.8185 

1.18 

3.0972 

4.2932 


.40 

2.7697 

3.3926 

+.9119 

1.20 

3.0998 

4.3148 


.42 

2.7975 

3.4192 


1.22 

3.1024 

4.3363 1 


.44 

2.8221 

3.4454 


1.24 

3.1049 

4.3579 


.46 

2.8442 

3.4712 


1.26 

3.1073 

4.3794 


.48 

2.8640 

3.4968 

+1.3500 

* 

1.28 

3.1097 

4.4009 


.50 

2.8818 

3.5220 

+1.4761 

1.30 

3.1121 

4.4224 


.52 

2.8980 

3.5469 

+1.6089 

1.32 

3.1144 

4.4439 


.54 

2.9127 

3.5715 


11.34 

3.1167 

4.4653 


.56 

2.9261 

3.5959 


1.36 

3.1190 

4.4868 


.58 

2.9384 

3.6201 

• 

1.38 

3.1212 

4.5082 


.60 

2.9497 

3.6441 


1.40 

3.1234 

4.5297 


.62 

2.9600 

3.6679 


|1.42 

3.1256 

4.5511 


.64 

2.9697 

3.6915 

+2.5419 

1.44 

3.1277 

4.5725 


. 66 

2.9786 

3.7149 

+2.7190 

1.46 

3.1299 

4.5939 


.66 

2.9869 

3.7382 

+2.9014 

1.48 

3.1320 

4.6153 


.70 

2.9946 

3.7614 

+3.0895 

1.50 

3.1341 

4.6367 


.72 

3.0019 

3.7844 

+3.2823 

1.52 

3.1363 

4.6581 


.74 

3.0087 

3.8073 


1.54 

3.1384 

4.6794 


.76 

3.0150 

3.8302 


1.56 

3,1405 

4.7008 


.78 

3.0210 

3.8529 


1.5706 

l 3.1416 

4.7124 

+10.177 
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Successive approximations to the function t(u) in the 
case of a circle (M^ = .406, A = .045) 


r 

CO 

f 0 (w) 

f X (^) 

f 2 (4») 


r 4 M 

f 5 (") 

de£. 

rad. 

u 

“TOOOOo 

.00000 

.00066 

"TOOOOO"' 

“7660 oo~ 

.00600 1 

.60000 i 

10 

.17453 

.17453 

.19144 

.18171 

.18611 

.19420 

.18558 

20 

.34907 

.34907 

.38086 

.36239 

.37082 

.36808 

.36976 

30 

.52360 

.52360 

. 56643 

.54122 

.55287 

.54900 

. 55135 

40 

.69813 

.69813 

.74683 

.71769 

.73138 

.72673 

.72951 

50 

.87266 

.87266 

.92136 

.89172 

.90587 

.90095 

.90388 

60 

1.04720 

1.04720 

1.09003 

1.06355 

1.07635 

1.07183 

1.07450 

70 

1.22173 

1.22173 

1.25352 

1.23362 

1.24333 

1.23985 

1.24191 

80 

1.39626 

1.39626 

1.41317 

1.40250 

1.40768 

1.40586 

1.40698 

90 

1.57030 

1.57080 

1 . 57080 

1.57080 

1.57079 

1.57080 

1.57080 
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Successive approximations to the function f (to) in the 
case of a symmetrical Joukowski profile (e = .15, 

M aP = .685, \= .469) 


W 




r 2 («) 


f 4 M 

jfljBJU 

racTT* 

0 

.0000 

.6600 

.OO00 

.0000“ 

rooocr 

.0000 

10 

.1745 

.0437 

.0324 

.0323 

.0324 

.0324 

20 

.3490 

.1273 

.1275 

.1274 

.1271 

.1275 

30 

.5236 

.2570 

.2786 

.2785 

.2779 

.2787 

40 

.6981 

.4470 

.4762 

.4755 

,4750 

.4760 

50 

.8726 

.6645 

.7090 

.7070 

.7068 

.7078 

60 

1.0472 

.9078 

.9649 

.9609 

.9615 

.9620 

70 

1.2217 

* 1.1682 

1.2320 

1.2257 

1.2270 

1.2272 

80 

1.3962 

1.4332 

1.4996 

1.4909 

1.4931 

1.4928 

90 

1.5708 

1.6956 

1.7586 

1.7484 

1.7514 

1.7508 

100 

1.7453 

1.9502 

2.0022 

1.9914 

1.9951 

1.9943 

110 

1.9198 

2.1857 

2.2253 

2.2156 

2.2193 

2.2182 

120 

2.0944 

2.4016 

2.4248 

2.4174 

2.4208 

2.4199 

130 

2.2689 

2.5922 

2.5992 

2.5956 

2 . 5982 

2.5972 

140 

2.4434 

2.7555 

2.7481 

2.7494 

2.7501 

2.7495 

150 

2.6180 

2.8902 

2.8728 

2.8787 

2.8776 

2.8778 

160 

2.7925 

2.9964 

2.9764 

2.9850 

2.9821 

2.9831 

170 

2.9670 

3.0773 

3.0639 

3.0702 

3.0675 

3.0686 

180 

3.1416 

3.1416 

3.1416 

3.1416 

3.1416 

3.1416 

_ 
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Velocity distribution about a circle (M^.406) 



Present 

Method 

Karman-Tsien 

Method 

— — ' 1 

Incompressible 

0 

0.000 

0.000 

0.000 

10 

.335 

.333 

.347 

20 

.675 

.667 

.684 

30 

1.014 

1.000 

1.000 

40 

1.350 

1.326 

1.286 

50 

1.671 

1.636 

1.532 

60 

1.952 

1.912 

1.732 

70 

2.185 

2.134 

1.879 

80 

2.336 

2.279 

1.970 

90 

2.389 

2.329 

2.000 


Table Illb 

Velocity distribution about a Joukowski Profile (£=.15, M^= 



Present 

Method 

Kdrm&n-Tsien 

Method 

Kaplan 

Incompressible 

0 


.849 

.839 

.870 

10 

.835 

.854 

.852 

.874 

20 

.847 

.869 

.867 

.887 

30 

.873 

.894 

.892 

.909 

40 

.912 

.928 

.927 

.938 

50 

.957 

.970 

.970 

.974 

60 

1.011 

1.019 

1.021 

1.016 

70 

1.073 

1.073 

1.078 

1.061 

80 

1.142 

1.132 

1.139 

1.109 

90 

1.215 

1.191 

1.203 

1.157 

100 

1.289 

1.250 

1.265 

1.203 

110 

1.360 

1.303 

1.322 

1.244 

120 

1.417 

1.347 

1.368 

1.278 

130 

1.446 

1.373 

1.394 

1.297 

140 

1.427 

1.369 

1.385 

1.294 

• 

150 

1.325 

1.308 

1.312 

1.247 

160 

1.088 

1.128 

1.110 

1 . 106 

170 

.645 

.704 

.674 

.738 

180 

0.000 

0.000 

0.000 

0.000 


.685) 




IACA TI to. 1006 


rig. 1 



Figure 1.- Velocity distribution along a circle (Moo*. 406). 
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